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Abstract. With the aim of clarifying the nature of the supernova events from stars having 
initial (at the main sequence) masses larger than ~ 8-10 Mq, we have developed a specif- 
ically tailored radiation hydrodynamics Lagrangian code, that enables us to simulate the 
evolution of the main observables (light curve, evolution of photospheric velocity and tem- 
perature) in these events. The code features and some test-case simulations as well as the 
possible applications of the code are briefly discussed. 
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1. Introduction 

It is widely accepted that core-collapse super- 
novae (CC-SNe) represent the final explosive 
evolutionary phase of stars having initial 
(i.e. at main sequ ence) masses larger than 
~ 8 -1 Mr. (see e.g.lHamuyl2003HHeger et all 
I2QQ3I: IWooslev & Weavej Il986l) . In addition 



to their intrinsic interest, these explosive 
events are relevant to many astrophysical 
issues associated, for example, with the nu- 
cleosynthesis processes of intermediate and 
trans-iron elements, the physical and chemical 
evolution of the galaxies, the production of 
neutrinos , cosmic rays a nd g ravitation al waves 
ICap pellaro & Turatto I 
'2003; Madau et al.' 

1989). 



(see e.g . lArnettI Il996t 



200C; 



1998 



iHaungs et al 
ilatnatunga & van den Bergh 



Moreover, CC-SNe seem to be particularly 



promising to measure cosmological distances, 
in addition to the well k nown me thod using 
type la SNe (see e.g. Nugent et al.l[2006, and 
references therein). 

Despite the importance of these explosive 
events in astrophysics, there are still ambigu- 
ities and uncertainties, linked to the extreme 
variety of CC-SNe displays (various energet- 
ics and different amounts of ejected material, 
leading to heterogeneous light curves, spec- 
tra and evolution of photospheric velocity and 
temperature) and related to the uncertainties 
in the modelling of stellar evolution and ex- 
plosion mechanism (see e.g. Filippenkoll 1 997t 
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iPumo et alJl2009HTurattol2063r and references 
therein). In particular the exact link between 
the physical properties of the SN explosion 
(ejected mass, explosion energy, stellar struc- 
ture and composition at the explosion) and the 
observational features is far from being well 
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Fig. 1. Effects of varying the ^*Ni mass on the light curve for models having the same total mass (16 Mq), 
initial radius (3 x lO'^cm) and total (kinetic plus thermal) energy (1 FOE = 10" ergs). The tracks correspond 
to 5'Ni masses equal to 0.005 Mq (red line), 0.010 Mq (green line), 0.020 Mg (blue line), 0.045 Mq (purple 
line), 0.055 Mq (cyan line), and 0.070 Mq (yellow line). 



established, and a "self-consistent" description 
of CC-SN events (from the quiescent evolu- 
tionary phases up the post-explosion evolution) 
is still missing. 

In this context, the creation of specific 
modelling tools that link progenitor evolution 
and explosion models to the main observables 
(i.e. light curve, evolution of photospheric ve- 
locity and temperature) of CC-SNe is of pri- 
mary importance for improving our knowl- 
edge about the nature of the CC-SN events, 
and the development of the new, radiation- 
hydrodynamics code, described in the next sec- 
tion, represents a key step in this direction. 



2. Code description and simulations 

The code, which solves the equations of 
relativistic, radiation hydrodynamics for a 
self-gravitating matter fluid interacting with 



radiation, is a new a nd improved ve r sion o f 
the one repor te d in Zampieri et al. (1996*), 
IZampieri etall (119981) and .Balberg et all 
(l2000h . 

Major strengt hs of the new vers ion of the 
code are (see also iPumo et al. I I2010I and Pumo 
& Zampieri 2010, in prep): 

- an accurate treatment of radiative transfer 
in all regimes (from completely optically 
thick to optically thin); 

- the coupling of the radiation moment equa- 
tions with the equations of relativistic 
hydrodynamics, adopting a fully implicit 
Lagrangian finite difference scheme, dur- 
ing all the post-explosive phases; 

- the computation of the evolution of the 
ejecta and the emitted luminosity taking 
into account the gravitational effects due to 
the compact remnant, enabling the code to 
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Fig. 2. Same as Figure[Tl but for the evolution of the photospheric velocity. 



deal with the fallback of material onto the 
compact remnant and, consequently, avoid- 
ing to overestimate the ejection of ^^Ni. 



These distinctive features make it possible 
to follow the post-explosion evolution of a CC- 
SN event in its "entirety", from the breakout of 
the shock wave at the stellar surface up to the 
so-called nebular stage, in which the ejected 
envelope has recombined and the energy bud- 
get is dominated by the radioactive decays of 
the heavy elements synthesized in the explo- 
sion. 

Consequently, the code is able to simulate 
the post-explosion evolution of the main ob- 
servables in CC-SN events for a given set of 
initial parameters (kinetic and thermal energy, 
envelope mass, radius, amount and distribution 
of ^^Ni) that are linked to the physical condi- 
tions of the CC-SN progenitor at explosion. 



3. Future developments 

We are starting to compute fine grids of mod- 
els evolved from post-explosion configurations 
and to carry out detailed investigations of the 
evolution of CC-SNe, in order to determine 
how the physical conditions of the CC-SN pro- 
genitor at the explosion affect the evolution of 
the main observables (see Fig.s[T]and|2]for the 
preliminary results of some of these calcula- 
tions). 

Our long term goal is the development 
of a sort of "CC-SNe Laboratory" interfacing 
our new relativistic, radiation-hydrodynamics 
code with stellar evolution, nucleosynthesis 
and spectral synthesis codes. This will allow 
us to describe the evolution of a CC-SN event 
in a "self-consistent" way from the evolution- 
ary stages preceding the main sequence up to 
the post-explosive phases as a function of ini- 
tial mass, metallicity and mass loss history of 
the CC-SN progenitor. 
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